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s»/ ■ We consider the two sequences of biorthogonal polynomials {pk,n)kLo 

r~| ! and {qk,n)kLo related to the Hermitian two- matrix model with poten- 
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of the coefficients in the recurrence relation satisfied by these polyno- 
mials, we obtain the limiting distribution of the zeros of the polyno- 
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1 Introduction 

1.1 Two- matrix model and biorthogonal polynomials 

The Hermitian two-matrix model is a probability measure 

J_g-nTr{y{Afi)+W(M2)-rMiM2) ^j^^ ^j^^ j--,^ j^ 

defined on the space of pairs of n x n Hermitian matrices (Mi, M2). Here, 
dMi dM2 is the standard Lebesgue measure on pairs of Hermitian matrices, 
V and W are the two potentials in the model which are typically polynomials 
of even degree with a positive leading coefficient, r / is a coupling constant 
and Zn is the normalizing constant 



With the two-matrix model (jl.ip . we associate two sequences of monic 
polynomials {pk,n)T=o and {qk,n)T=Q, where degpk^n = degqk,n = k, called 
biorthogonal polynomials, defined by the property 

Pfc,n(x)g,,„(y)e-"(^(^)+^(^)--^^) dx dy = 0, k^ j. (1.2) 

Existence and uniqueness of these polynomials was proved by Ercolani and 
McLaughlin in [9]. They also showed that the zeros of these polynomials 
are real and simple. It will be one of our results that the zeros interlace, see 
Theorem 12.11 

It is well-known that the eigenvalue correlations of the matrices Mi and 
M2 from p.ip are determinantal with correlation kernels that can be ex- 
pressed in terms of the biorthogonal polynomials and their transforms, see 
[21 [3l El [ini E] . As an example of these relations we have that 

E[det{xln - Ml)] = pn,n{x), 

E[det{yln-M2)]=qn,n{y)^ 

which show that the diagonal biorthogonal polynomials pn,n and g„^„ can 
be considered as 'typical' characteristic polynomials for Mi and M2, respec- 
tively. In this sense the zeros of Pn^n are typical eigenvalues of the matrix 
Ml. This explains our interest in these zeros as n — )• 00. 



1.2 A vector equilibrium problem 

In [S] two of us studied the limiting eigenvalue behavior of the matrix Mi 
in the two-matrix model (jl.ip for the case of an even polynomial V with 
positive leading coefficient and 

W{y) = \y\ 

The biorthogonal polynomial "Pn^n associated with this model can be charac- 
terized by a Riemann-Hilbert problem of size 4x4 [3l[l2]. The Deift-Zhou 
steepest descent method was successfully applied to the Riemann-Hilbert 
problem from |12j . A crucial ingredient in [8] is the introduction of a vector 
equilibrium problem with external field and upper constraint that describes 
the limiting mean eigenvalue distribution. 

To state the vector equilibrium problem we use notions from logarithmic 
potential theory |16] . We define the logarithmic energy of a finite positive 
measure i/ on C as 

^(^) = //logl -Av(x)^v{y), 

JJ \x-y\ 

and the logarithmic potential as 

U^iz) = /log , ^ , dt/(x). (1.3) 

J \z-x\ 

If 1^1 and 1^2 are positive measures on C with I(z^i), 1(1^2) < 00, we also define 
their mutual logarithmic energy as 

1(1^1,1^2) = / / log -. r di/i(x) dz^2(y)- 

JJ \x-y\ 

The vector equilibrium problem from [8] asks to minimize the energy 
functional 

/(pi) - I{pi,p2) + I{P2) - Iip2,P3) + HPs) + I (V{x) - ^\tx\'/A dpiix), 

among all measures pi , p2 and p^ with finite logarithmic energy that satisfy 

(a) pi is supported on M and /5i(M) = 1; 

(b) p2 is supported on iR. and p2{iR) = 2/3; 

(c) p3 is supported on M and /33(M) = 1/3; 



(d) p2 satisfies the constraint P2 < o" where a is the unbounded measure 
on zM defined as 

27ri 
Here, dz is the complex hne element on iM. 

In [8] it is shown that there is a unique minimizer {I'l, 1^2,1^3) of the 
vector equilibrium problem. The measure vi is supported on a finite union 
of disjoint intervals. For the case of one interval (one-cut case) it was shown 
in [8] that the density of z^i is the limiting mean density of the eigenvalues 
of Ml. See [15] for the extension to the multi-cut case. 



1.3 Aim of this paper 

It is the aim of this paper to give a new perspective on the nature of the 
above vector equilibrium problem. In [8] the vector equilibrium problem 
was simply posed out of the blue, while in the present work it arises after 
certain calculations. 

In addition, we derive a similar vector equilibrium problem for the case 
that 

Vix) = ^ and P7(y) = ^+t^, (1.4) 

where i G M is a real parameter. When t < this vector equilibrium problem 
has the novel feature that an external field is acting on the third measure 
as well. 

Due to the fact that V{x) is quadratic in (|1.4p the second sequence {qk,n)k 
of biorthogonal polynomials is actually a sequence of orthogonal polynomials 
that satisfy a three term recurrence relation. The biorthogonal polynomials 
Pk^n also satisfy a recurrence relation with recurrence coefficients that can 
be expressed in terms of the recurrence coefficients for qk,n- We are going 
to analyze these recurrence coefficients as n — )• 00 and obtain the vector 
equilibrium problem from this analysis. This approach is similar in spirit 
to the analysis of [13j for the recurrence coefficients of multiple orthogonal 
polynomials in a model of non-intersecting squared Bessel paths. 

Furthermore, we prove that the first component z^i of the vector of mea- 
sures (z^i, 1^2; 1^3) minimizing this vector equilibrium problem is equal to the 
limiting zero distribution of the diagonal polynomials pn,n as n — )• cx). The 
measure vi is also the limiting mean eigenvalue distribution of the matrices 
Ml in the two-matrix model. This will be proved in a forthcoming paper. 

In the next section we give a more detailed description of the main results 
in this paper. 



2 Statement of main results 

2.1 Interlacing zeros 

Our first result deals with biorthogonal polynomials defined by (jl.2p with 
general potentials V and W. Our result is that the zeros of consecutive 
biorthogonal pfc,n and Pk+i,n are interlacing. This was proved by Woerdeman 
[IB] for a special case. Two ordered sequences of real numbers ai, . . . ,ak 
and /?!,..., /3fe+i are said to interlace if 

(3i <ai < 132 <a2 < ■■■ < au^i < Pk < ak < A+i- 

Theorem 2.1. Take r 7^ and suppose that V and W are functions for 
which the integrals 

/OO /'OO 
/ xS^e-'^(^(^)+^(^)-^^^)dxd2/, A; = 0,1, 2,..., i = 0,1,2,..., 
-00 J —00 

converge. Define the two sequences ofmonic biorthogonal polynomials {pk,n)'kLo 
and iqk,n)kLo ^^ '^^ (|1.2p . Then the following statements hold for every 
k = l,2,... 

(a) The zeros of pk,n o-nd pk+i,n interlace, and similarly, the zeros of qt^n 
and qk+i,n interlace. 

(b) If the potentials V and W are even, then the positive zeros ofpk^n o-nd 
Pk+2,n interlace, and similarly, the positive zeros of qk,n CLnd qk+2,n 
interlace. 

Theorem 12.11 is proved in Section [3j 

2.2 Limit of zero counting measures 

In the rest of the paper we restrict ourselves to the quadratic and quartic 
potentials (II. 4p . The biorthogonal polynomials associated with this model 
are thus defined by 



00 /•OO 



00 J —00 



Pfc,n(x)<Zi,n(?/)e-"("'/2+^'/^+*^'/2-™^) dx dy = 0, k^ j. 

(2.1) 



We also assume without loss of generality that 

r >0. 



As a second result, we will show that in this case the limiting zero dis- 
tribution of the diagonal polynomials pn,n exists. For a polynomial P of 
degree n, we introduce the normalized zero counting measure 



K^) = ^ E '^^ 



n 

P{x)=0 



where the sum is taken over all zeros of P counted with multiplicity. We 
say that a sequence of measures m„ converges weakly to the measure v if 

fdun^ If dz/, 



for every bounded continuous function /. 

Theorem 2.2. There exists a Borel probability measure vi with supp(i^i) C 
R such that ui is the limiting zero distribution of the diagonal polynomials 

lim u{pn,n) = J^l, 
n— J-oo 

where the limit is in the sense of weak convergence of measures. 

The proof of Theorem 12.21 is given in Section [71 It is constructive. For 
an explicit formula of the limiting zero distribution i^i , we refer to (jT.lh . 

Note that Theorem 12.21 is about the limiting distribution of the zeros of 
the biorthogonal polynomial Pn,n- The measure z^i is also the limiting mean 
distribution of the eigenvalues of the matrix Mi from the two-matrix model, 
but this is not proved in this paper. This will follow from an analysis of the 
Riemann-Hilbert problem as in [8J, which is under current investigation. 

2.3 Vector equilibrium problem 

Our final result is that the limiting zero distribution z^i can be characterized 
by a vector equilibrium problem depending on external fields Vi and V3, 
and on a constraint a. These objects will be described next in terms of the 
solutions of the equation 

u;^ + tuj = Tz. (2.2) 

For t = 0, the vector equilibrium problem reduces to the vector equilibrium 
problem described in Section 11.21 



External field Vi on M. For z = x gM., the equation ()2.2p has either one 
or three real solutions. We use u}i{x) to denote the real solution with the 
largest absolute value. This is the only real solution for t > and also for 
t < and z = X with |3;| > x* where 

2f-t)3/2 
X* = x*(t) = ^-^ — , t<0. (2.3) 

3\/3r 



For t < and —x* < x < x* there are three real solutions of (j2.2p which we 
denote by ujj(x), j = 1, 2, 3, and which we number such that 

|a;i(2;)| > \u2{x)\ > \i03{x)\. (2.4) 

This means in fact that 0J2{x) < uj^^x) < < uJi{x) if x G (0,x*) and 
iOi{x) < < ujsix) < uj2{x) if X G (—X*, 0). 

In both cases, the external field Vi is defined by 

x2 
Vi{x) = — + min(VF(y) — rxy) 



= y--a;i(x)4--ta;i(x)2, x G R. (2.5) 

The second identity in ()2.5p comes from the fact that the minimum of W{y) — 
Txy is taken at y = iOi{x) and rx = a;i(x)^ + tuji{x) by (12. 2p . 

External field V3 on M The external field V3 vanishes identically for t > 

Vsix) = 0, for X G M, if t > 0. (2.6) 

For t < and x G M, we define 

y(^) ^ /i^2(x)^ + ita;2(a;)2-|a;3(x)4-itw3(x)2, for |x| < x*, 
1 0, for |x| > X*. 

where x* is given in (|2.3p . and uj2{x) and 0^3 (x) are the solutions of (|2.2p 
that satisfy (ITil) . 

While Vi(x) is related to the global minimum of the function 

W{y) -Txy = ^ + t^- Txy, y G M, (2.8) 

V3(x) can be interpreted as the positive difference between the local maxi- 
mum and the other local minimum of (j2.8p on M, which indeed exist if and 
only if f < and |x| < x*. 



Upper constraint a on iM For z = iy £ iM, the equation (j2.2p has either 
one or three purely imaginary solutions. There are three purely imaginary 
solutions if and only if t > and |y| < y* where 

2*3/2 
y* = y*(t) = ^-^, t>0. (2.9) 

Otherwise there is only one purely imaginary solution and the two other 
solutions are located symmetrically with respect to the imaginary axis. We 
then let uji{z) be the solution of (j2.2p with positive real part. For conve- 
nience we put 

y* = 0, if t < 0. 

The upper constraint a is defined as follows. The support of a is 

supp(cr) = (— ioo, —iy*] U [iy* , ioo), (2.10) 

and a has the density on supp(o") given by 

— Y = -Rewi(z), zesupp(a), (2.11) 

1 0.2 1 vr 

for every fixed t G M. 

Now we state our final main result. 

Theorem 2.3. The measure ui from Theorem \2.2\ is the first component of 
the unique vector of measures {1^1,1^2,1^3) minimizing the energy functional 

E{pi,P2,P3) =I{pl) -I{PliP2)+I{P2) - I{P2,P3) + I{P3) 

+ fvi{x)dpi{x)+ fv3{x)dps{x). 



among all vectors {pi, P2, P3) of measures with finite logarithmic energy sat- 
isfying 

(a) pi is supported on M and pi{M.) = 1, 

(b) p2 is supported on iM and p2{i^) = 2/3, 

(c) p3 is supported on M and /53(M) = 1/3, 

(d) p2 satisfies the constraint P2 ^ <y- 



Here Vi and V3 are defined in (j2.5p and (|2.6p - (|2.7p . respectively, and a is 
defined by (121^1) - (l2TT]l . 

Theorem 12.31 is proved in Section [HI 



2.4 Phase diagram and critical behavior 

The proof of Theorem [213] is constructive, and we find fairly exphcit formulas 
for the minimizing measures Uj. Indeed, we obtain i/j as an average 



of measures fi- depending on a parameter ^ and these measures are given 
by formulas (j5.15p and (|6.15p . 

It follows from the analysis leading to these formulas that the supports 
of the measures vi, a — 1^2, and v^ have the following form 

supp(z^i) = [-a, -13] U [13, a], 
supp(cr - 1^2) = iK \ (-^7, i-y), 
snM^3)=R\ (-6,6), 

for some a>/3>0, 7, 5>0 depending on t S M and r > 0. 

We may distinguish a number of cases, depending on whether (3, 7, or 5 
are equal to zero, or not. At least one of these is zero, and generically, no 
two consecutive ones are zero. Our analysis leads to the phase diagram in 
the tr-plane shown in Figure [H 

Case I: /5 = 0, 7 > 0, and 6 = 0. Thus in this case there are no gaps in the 
supports of the measures 1^1 and 1^3 on the real line. The constraint a 
is active along an interval [—17, ij] on the imaginary axis. 

Case II: /3 = 0, 7 > 0, and 6 > 0. In this case the measure vi is still 
supported on one interval. However there is a gap {—S,6) in the sup- 
port of f3. As in Case I, the constraint a is active along an interval 
[— ^7, 47] on the imaginary axis. 

Case III: /3 > 0, 7 = 0, and 6 > 0. In Case III there is a gap in the 
supports of vi and ^3, but the constraint on the imaginary axis is not 
active. 

Case IV: /3 > 0, 7 > 0, and 6 = 0. In Case IV there is a gap in the 
support of vi, but there is no gap in the support of z/3, which is now 
the full real line. The constraint is active along an interval along the 
imaginary axis. 

Cases II and III are new in the sense that they do not appear in [8] . The 
opening of a gap in the support of v^ is due to the external field V3 that acts 



on the third measure in the vector equihbrium problem. As V3 is identically 
zero for t > 0, the Cases II and III do not appear if t > 0, as can be seen in 
Figure [H 

Critical behavior occurs at the curves that separate the different cases 
from each other. These critical curves are given by the equations 

T = Vt + 2, -2 < t < 00, and '^ = \ > -00 < t < 0. 

On the critical curves two of the numbers /3, 7 and 5 are equal to zero. For 
example, on the curve between Case II and Case III, we have /3 = 7 = 0, 
while 6 > 0. Finally, note the multi-critical point 

t = -l, T = l 

in the phase diagram, where /3 = 7 = (5 = 0. All four cases come together 
at this point in the tr-plane. 

We do not discuss the critical and multi-critical behavior any further 
in this paper. However, it would be particularly interesting to analyze the 
nature of the multi-critical point. 

2.5 Overview of the rest of the paper 

The rest of the paper is devoted to the proofs of the three main results. 
Theorem 12. II is proved in Section [3l The other sections deal with the specific 
model ()1.4p and lead to the proofs of Theorems 12. 21 and 12. 31 in Sections [7] and 
IHl The intermediate sections contain auxiliary results that will be essential 
to these proofs. 

In Section 0] it is shown that the biorthogonal polynomials satisfy re- 
currence relations with recurrence coefficients that have certain asymptotic 
behaviors. We introduce a new parameter ^ and consider the asymptotic 
behavior in the t^-phase space. Two different types of asymptotic behavior 
will lead to a separation of the phase space into two regions Ci and C2, 
where Ci is the one-cut case region, and C2 the two-cut case region. 

With every ^-value, we associate a vector equilibrium problem for three 
measures. These equilibrium problems serve as building blocks for the vector 
equilibrium problem of Theorem 12.31 We introduce and analyze the one-cut 
case in Section [5] and the two-cut case in Section [6l Due to the different 
asymptotic behavior of recurrence coefficients, the analysis in both sections 
is significantly different. 
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Figure 1: The phase diagram in the tr-plane: the critical curves r = \/T+^ 
and T = \/ —J separate the four cases. The cases are distinghuished by the 
fact whether is in the support of the measures vi, o — v^^ and 1/3, or not. 

3 Proof of Theorem [23] 

In this section we prove Theorem 12.11 The proof of this theorem is inspired 
by [1] and uses the following theorem from [9]. 

Theorem 3.1. Suppose that r 7^ 0. Let V and W be functions for which 
the integrals 

fOO />oo 

3.fcyie-n(v(x)+H/(y)-rxj,)^^^y^ A; = 0,1,2,..., j = 0,1,2,..., 



—00 J —00 



converge. Then, for each k, there is a unique monic polynomial pk,n o.nd a 
unique monic polynomial Qk^n such that the families of polynomials {pk,n)'kLo 
and {qk,n)'kLo satisfy ( tJ.^j) . In addition, the zeros of these polynomials are 
real and simple. 

Proof This is P Theorem 1]. D 

To establish Theorem l2.1l fa) it is clearly enough to prove the statements 
about the zeros of Pk,n, since the results about the zeros of qk,n follow by 
symmetry. 
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Proof of Theorem \2.1\ (a). Fix an integer k>\ and consider the linear com- 
bination Apk^n~\~Bpk+i^n with [A, B) 7^ (0, 0). We claim that this polynomial 
has no real multiple zeros. 

To see this, assume that xq is a real zero of multiplicity at least two. 
Then we can write 

Apk,n{^) + Bpk+iA^) = {x- XQfr{x), 

where r is polynomial of degree < k — 1. From the biorthogonality (|l.'2p it 
follows that 

r{x)y\x - 2;o)2e-«(^W+W^fe)-™3/) dxdy 

r(2;)y'e-"(^(^)-i^°§l^-^ol+^(j')-™3')dxd2/ = 0, 



for / = 0, 1,2,...,A; — 1. Applying the uniqueness part of Theorem 13.11 
to the modified potential V{x) — -loglx — xo|, we obtain r = 0. Thus, 
A = B = Q^ which yields a contradiction and, therefore, proves the claim 
that Apk^n + Bpk+i,n has no real multiple zeros if (A, B) / (0, 0). It follows 
that the linear system of equations 

'Pk,n{^) Pk+l,n{x)\ fA\ ^ /O 
P'k,rX^) P'k+l,ni^)J \Bj {0 

has only the trivial solution A = B = 0, for every x G M. Therefore the 
matrix has a non-zero determinant, and thus 

Pk,n{x)p'k+l,nix) - p'k,n{x)pk+l,n{x) 7^ 0' X G M. 

By continuity and the behavior as x — t- oo, we conclude from this that 

Pk,n{^)p'k+l,n{x) - p'k,n(.^)Pk+l,n{^) > 0' X G M. (3.1) 

Now consider two consecutive zeros xi and x/+i of p^+i.n- Because these 
zeros are simple, we have that 

P'k+l,ni^l)p'k+l,ni^l+l) < 0- 

From ([XI]) we find Pk,n{^l)p'k+i^n(^l) > and pk,n{xi+i)p',^^^^^{xi+i) > 0. 
Hence, we obtain 

Pk,ni^l)Pk,n{^l+l) < 0. 

Therefore, pk^n must have a zero between xi and xi^i. Hence, in between 
any two consecutive zeros of Pk+i,n, there is a zero of Pk,m which implies 
that the zeros of pk,n and Pk+i,n interlace. D 

12 



Proof of Theorem \2.1\ (h). The proof of (b) follows the same strategy. Let 
A; > 1 and assume that the linear combination ^Pfc,n + ^Pfc+2,n has a positive 
multiple root, say xq > 0. Because V and W are even potentials, pk^n and 
Pk+2,n are either both even or both odd. Therefore, also —xq is a double 
zero and we can write 

where r is polynomial of degree < k — 2. Then, r = and, thus, A = B = 
as in the proof of part (a). 

Hence Ap^^n + ^Pfc+2,n with (^4, B) ^ (0, 0) has no positive multiple 
zeros. Therefore, the linear system of equations 

Pk,n{^) Pk+2,n{^)\ I'A'^. _ /O 
P'kJ^) P'k+2,ni^)) \BJ - \0 

has only the trivial solution for every x > 0. Thus, as in the proof of part 

(a), 

Pk,n{x)p'k+2,nix) - Pk^n{x)Pk+2,n{x) > '^^ X > 0. (3.2) 

The proof of interlacing of the positive zeros follows from p.2p in the same 
way as before. D 

Remark. It was shown in [9] that the zeros of biorthogonal polynomials 
are real and simple. We can prove this result in an alternative way as follows. 
First assume that pk^n has a non-real zero xq = a -\- bi, a,b £ M, b ^ 0. 
Then also xq = a — bi is a zero of Pk,n- Thus, pk^n can be written as 

Pk,n{^) = {{x - af + b^) r(x), 

where r is a polynomial of degree k — 2. Now observe that 

y' {{x - af + b^) r(x)e-"(^(^)+^(j^)-™^) dxdy 

(x)y'e-"(^(")-^^°g(("-")'+''')+^(^)-™J') dxdy = 0, 



oo 

oo roo 



r[ 

oo 



for I = 0,1,2, . . . ,k — 1. Note that the second equality follows from (jl.2p . 
Applying the uniqueness part of Theorem 13. II to W{y) and the modified po- 
tential V{x) log ((x — af + 6^) we obtain r = which is a contradiction. 

Therefore the zeros of pk^n are real. Moreover, by putting 6 = in the above 
argument we also obtain that the zeros are simple. 

Of course, a similar reasoning shows that the zeros of qk,n are real and 
simple. 
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4 Preliminaries on recurrence coefficients 

4.1 Relations between recurrence coefficients 

In the rest of the paper we consider the model with quartic and quadratic 
potentials (jl.4p . In that case both sequences {pk,n)k and {qk,n)k of biorthog- 
onal polynomials defined by (|2.ip satisfy a recurrence relation with recur- 
rence coefficients that are related to each other as described in the following 
lemma. 

Lemma 4.1. For n > 1, the sequence of polynomials {qk,n)'kLo ^s orthogonal 
with respect to the weight 

nj{y)=e-<y'/^+'y'/^--'y'/^\ yeW, (4.1) 

and, therefore, satisfies a recurrence relation of the form, 

yQk,n{y) = Qk+iAv) + ak,nqk-i,n{y)^ A: = 0, 1, 2, ... , 

where q~i^n = 0, ao,n = 0, and au^n > 0, /c = 1, 2, . . . . 

In addition, the sequence of polynomials iPk,n)'kLo satisfies a recurrence 
relation of the form 

XPk,n{^) =Pk+l,n(x) + ^k,nPk-l,n{^) + Ck,nPk-3,n{^), k = 0,1,2,... , (4.2) 

where p_i,n = P-2,n = P-3,n = 0. The recurrence coefficients are related as 
follows 

k 

h,n = ak,n{ak~l,n + ak,n + ak+l,n + t) =T ak,n+ -, k>l, (4.3) 

n 



Ck,' 



T^afc-2,n«fc-i,nafc,n, A; > 3. (4.4) 



Proof. The lemma is known in much greater generality for general polyno- 
mial potentials V and W, see [2]. 

Explicit calculations that lead to (|4.2p . (j4.3p and (14. 4[) for the case t = 
are given in |18j . These calculations extend to general t G M in a straight- 
forward way. D 

4.2 Asymptotic behavior of recurrence coefficients 

Our next result concerns the asymptotic behavior of the recurrence coef- 
ficients ak,n, bk,n and Ck,n as A:,n — )■ oo such that k/n — )• ^ > 0. Let us 
introduce some notation. We write 

lim Xk.n = X, 

k/n^>S, 
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if lim Xk- n- = X holds for every two sequences of positive integers (kj) 

and (rij) that satisfy kj,nj — t- c« and kj/uj — >■ ^ as j — >■ oo. In the same 
spirit we write 

hm Xk^n = X 

k even 



if hm X^. n- = X holds for every sequence of positive even integers (A;,) 

and every sequence of positive integers {uj) that satisfy kj^rij — )• cxd and 
kj/uj — )• ^ if j — )• oo. The limit with the subscript 'odd' is defined similarly. 
The limiting behavior of the recurrence coefficients a^^^, 6^^^, Cfc^„ as 
/c, n — )■ oo, k/n — >• ^ depends on the values of t G M and ^ > 0. We consider 
the coupling constant r > to be fixed. We define the critical ^-values 

In the t^-plane the equation C = Ccr, t < t^, defines a semi-parabola that 
separates the upper half of the t^-plane into two regions 

Ci ■■ S, > Ccr, 

C2 : < ^ < ^cr, -00 < t < T^, 

see Figure [2j We refer to Ci as the one-cut case region, since the zeros of 
the orthogonal polynomials qk^n accumulate on one interval as A;, n — >• 00 
and k/n — )■ ^ > .^cr- If t < r^ and ^ G (0,^cr) the zeros of qk,n accumulate 
on two disjoint intervals and therefore we call C2 the two-cut case. 
We now state the following theorem. 

Theorem 4.2. (a) //^ > ^cr then the limits of ak^n, ^k,n, Ck,n o,s k/n — t- ^ 
exist and we have 

hm ak,n = a(6 := ^^—^ , (4.5) 

k/n->( 

hm 6fc,„ = b{0 := aiOiMO +t)= r^a{0 + ^, (4.6) 

fc/n— >-^ 

hm Cfc_„ = c(e):=rV(e). (4.7) 

fc/n— >■$ 

(b) If t < t'^ and < ^ < ^cr; then the recurrence coefficients ak,n, bk,n, 
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Figure 2: t,^-phase diagram: the semi-parabola separates the one-cut case 
region Ci from the two-cut case region C2. 



Ck,n exhibit 2-periodic behavior as k/n — )■ ^ and we have 



r^-t-y^i^- 


-ty- 


-Ai 


2 


r'-t+y/ir^- 


-t?- 


-Ai 



lim ak,n = ao(0 - 

k even 

hm ak,n = ai{i) := 

k odd 

lim bk,n = boiO ■■= aoiOiaoiO + 2ai(0 + *), 

k/n—^t; 
k even 

lim bk,n = biiO ■■= ai(0(2ao(e) + ai{0 + t), 

k odd 

lim Ck,n = co(0 := T^al{^)ai{^), 

k even 

lim Ck,n = ci{i) := T^ao(Oai(0- 

k/n^t^ 
k odd 



(4.8) 

(4.9) 
(4.10) 
(4.11) 
(4.12) 
(4.13) 



(c) If t < t"^ and C = ■^cr; then 

a(e) = ao(0 = ai(0, HO = bo{0 = bi{0 , c(0 = co(0 = ci(0, 

(4.14) 

and all of the above limit relations continue to hold for ^ = .^cr- 

Proof. The recurrence coefficients ak,n appear in a recurrence relation for 
the orthogonal polynomials pk,n- These polynomials are orthogonal with 
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respect to the weight (j4.ip . For this type of orthogonal polynomials the 
asymptotic behavior of recurrence coefficients was studied by Bleher and Its 
in two papers. The paper |1] deals with the two-cut case < ^ < ^cr- In 
[5] the critical case ^ = ^cr is studied and the results for the one-cut case 
^ > .^cr are given as well. 

Having ([33]), gSD, and gj]), the limits, (gSD, (|32D, and ()i30]) - pj3D 
follow directly from Lemma l4. II D 

5 Asymptotic analysis in one-cut case 

5.1 Results from the literature 

In what follows we will associate with each ^ > a function of the form 

s{w) = w + dW + ^ + ^ + ^, d(3) ^ 0. (5.1) 

Such functions appear as symbols of banded Toeplitz matrices [6], and we 
need certain results [3 [TH] that were derived in that context. Although we 
will not use Toeplitz matrices in this paper, we still refer to s as the symbol. 
We denote the solutions of the algebraic equation s{'w) = z by Wj{z), 
J = 1, . . . , 4 and order them by their absolute value, such that 

\wi{z)\ > \W2(,Z)\ > \W3{Z)\ > \W4{Z)\ > 0. (5.2) 

Typically, there is strict inequality in ()5.2p . If for certain z G C two solutions 
have the same absolute value, then we pick an arbitrary numbering that 
satisfies (|5.2p . Furthermore, we define 

rj = {zeC\\w,{z)\ = \wj+i{z)\}, j = 1,2,3, (5.3) 

which are finite unions of analytic arcs and exceptional points, see [6l [7]. 

A point z E C for which the algebraic equation s{w) = z has a multiple 

solution is called a branch point. 

We use the solutions Wj{z) to the algebraic equation to define three Borel 

measures 

1 4^(w' {z) w'{z)\ 
d^,(z) = — y ^^ - ^±^4 dz, (5.4) 

^^^ > 2Tli ^^\wk-{z) Wk+{Z)J ^ ' 

for z G Tj, j = 1,2,3. Here, it is assumed that every analytic arc of Tj is 
provided with an orientation and that dz denotes the complex line element 
on Tj according to this orientation. Furthermore, Wk±{z) is the limiting 
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value of Wk{z) as z ^ z from the it side on each of the arcs in Tj. The 
+ side (— side) is on the left (right) if one traverses Tj according to the 
orientation. 

The vector of measures (/^i,/^2i/^3) is characterized as the unique mini- 
mizer of a vector equilibrium problem. 

Theorem 5.1. Define the energy functional Eq as 

Eo{Pl,P2,P3) =I{Pl) -I{Pl^P2) + I{P2) - I{P2,P3) + I(.P3), 

where pi,P2,P3 olt^ positive measures on C with finite logarithmic energy. 
Then the following statements hold. 

(a) The vector of measures (pi, P2, Ps) given by (j5.4p is the unique min- 
imizer for the functional Eq among all vectors {pi,P2-,Pz) of positive 
measures with finite logarithmic energy, satisfying 

(i) supp(pj) C Tj, for i = 1,2,3, and 
(ii) pi{Ti) = I, p2{T2) = 2/3, and palFs) = 1/3. 

(b) The measures pi, P2i P3 satisfy for some constant i 

wi{z) 



i - lUf"' (z) + V^z) =log 

U^^'iz) - 2Uf'^z) + U^'-'iz) = log 

U^^{z)-2U^''{z) =log 



W2{z) 
W2{z) 



wsiz) 

W3{z) 



■W4[z) 



(5.5) 
(5.6) 
(5.7) 



for every z £ C 



Proof. The proof of Theorem 15.11 can be found in [7]. The conditions 
(j5.7p are the Euler-Lagrange variational conditions for the vector equilibrium 
problem. Note that the right-hand side of the jth variational condition 
vanishes if z S Tj . In [7j there also appear constants £2 and £3 in ()5.6p and 
(|5.7p . However, these constants vanish because r2 and T3 are unbounded 
and 

[/^^■(z) = -//,(r,)iog|z| + o(i). 



as z 



00. 



D 



If the symbol ()5.ip depends on a parameter .^ > 0, say 



.„,<..«) + ;!lM + iq« + ^, 



w 



w^ 



w^ 



(5. 
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then we use Wj{z; ^), rj(^) and fi, to indicate the dependence of the notions 
from ()5.2p . (j5.3p and ()5.4p . respectively, on the parameter ^. 

Next, we state a result of Kuijlaars and Roman [13] on polynomials 
satisfying certain recurrence relations. It will be the key ingredient of the 
proof of Theorem 12.21 given in Section [71 

Theorem 5.2. Let for each n € N a sequence ofmonic polynomials {pk,n)'kLo 
be given where deg pk,n = k. Furthermore, suppose that 

(a) these polynomials satisfy the recurrence relations 

XpU-) = P.+l,n(x) + 4V.n + 4'1P.-I,n + 4V-2.n + 4V-3.n, 

(5.9) 

for certain real recurrence coefficients d\. !^, j = 0, 1, 2, 3; 



""k 



(b) the polynomials pk^n have real and simple zeros x^'^ < ■ ■ ■ < x^^'" 
satisfying for each k and n the interlacing property 

fc+l,n , k,n , k+l,n s ■ i i 

Xj <Xj' < Xj^{ , for J = l,...,k; 

(c) for each j = 0, 1, 2, 3 the set of recurrence coefficients 

{dgi|/c + l<n} 
is bounded; 



(d) there exist continuous functions d^^' : (0, +oo) — t- M, j = 0, 1, 2, 3, such 
that for each .^ > 

lim 4^'i = dO)(0, (5.10) 

andd^^\C) /O; 

(e) we have 

ri(^)cM, for every C>0 

where Ti{^) is the set defined as in (|5.3p corresponding to the ^- 
dependent function ([5.8p with d^^' (^) coming from ([5.10p . 

Then, the normalized zero counting measures i'{pk,n) have a weak limit 
as /c, n —>• oo with k/n — )■ A > given by 

1 r^ 

lim iy{pk,n) = ^ fJ'id^, (5.11) 

k/n^fX A Jq 

where for each ^ > the measure fi\ is given by ([5.4p corresponding to the 
function (|5.8p . 
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Proof. This is [13^ Theorem 1.2] for the case of a five term recurrence ()5.9 



D 

The intuition behind Theorem l5.2l is that the zeros oipt^n are eigenvalues 
of a A; X A; matrix with five non-zero diagonals 



' "0,n 



1 



"l,n "l,n 







\ 



(2) ,(1) ,(0) 



^2,n 
"3,n 





"2,n 
"■3,n 



^2,n 
"3,n 





1 

.(0) 
^3,n 





1 



V 



(5.12) 



• / 



Under the assumption (|5.10|) the entries are slowly varying along the di- 
agonals if k and n are large, so that locally the matrix (j5.12p looks like 
a five-diagonal Toeplitz matrix. Then for each ^ one considers the exact 
Toeplitz matrices with the entries d^'^\i) , (P'\^) , (P-\C) ■, d'"^\i) ■, 1 along the 
diagonals for which it is known, see [Si [71 Ull |T7] , that the eigenvalues accu- 
mulate on ri(^) as the size grows, with ijl\ as limiting normalized eigenvalue 
counting measure. The distribution of the eigenvalues of (j5.12p is then ob- 
tained by averaging of the measures ^\ as in (jS.lip . 



5.2 Analysis of the symbol si in the one-cut case 



It will be our goal to apply Theorem 15.21 to the biorthogonal polynomials 
Pk^n that have the recurrence relation (j4.2p . The recurrence coefficients in 
(j4.2p have the appropriate limits only in the one-cut case. We discuss this 
case first. 

We therefore assume that E, > Ccr- In that case we have by Theorem 
14.21 that the recurrence coefficients hk^n, Ck,n have limits 6(^) and c(^) given 
by (|4.5p and (|4.7p as A:, n — )• cx) and k/n — )• ^. We therefore associate with 
C > '^cr the symbol 



si{w;C) = w 



+ 



w 



w^ 



As already noted before, we use Wj{z;^) for j = 1,2,3,4, and Tj{^). 
j = 1, 2, 3 to denote the quantities related to the symbol (j5.13p . 



(5.13) 



^f for 
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In order to apply Theorem 15.21 we need to know that ri(^) C M, see 
assumption (e) in Theorem 15.21 Then the proof of Theorem 12.31 fohows the 
approach outhned in [131 Section 7]. So, to obtain an external field Vi acting 
on j/i and an upper constraint a acting on 1^2 we will need that ri(^) is an 
increasing and r2(^) a decreasing set as a function of ^. These features are 
contained in the following theorem. 

Theorem 5.3. Let r > and t G M. Then for every S, > £,ct we have that 

ri(^) C M, T2iO C iR, and r3(e) C R. 

For ^ > ^cr the set Ti{S^) is increasing as a function of ^, while the set T2{Cj 
is decreasing. More precisely, there exist a('C))7('C) > such that 

ri(e) = [-a(e),a(e)], 

r2(e) = iM\H7(e),n(6), (5.14) 

rm = R- 

In addition, we have that 

(a) C '"^ ck(C) ^•s strictly increasing for ^ > ^cr with 

lim a(^) = +00, and lim a(^) = if t > t , 

(b) ^ I— 7- 7(^) is strictly increasing for ^ > ^cr with 

lim 7(^) = +00, and lim 7(^) = y* if t > t , 

see (12. 9p for the definition of y* . 

Figure [3] shows the sets ri(,^), r2(^), and r3(^) in the complex plane in 
the one-cut case. 

The proof of Theorem 15.31 is given in the next subsection. Here we note 
that as a consequence of Theorem 15.31 two consecutive sets among ri(^), 
r2(^) and r3(^) are not overlapping, which implies that (j5.4p can be written 
more simply as 

. , 1 /K„(^;0 w'j (z;^)\ 



dfi^z) = — '- - '^ dz, z G r,(C), j = 1, 2, 3. 

^ 2m \Wj_{z;C) Wj^{z;^) I 



(5.15) 
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, r2(0 
ri(e) 1 



r3(0 



Figure 3: The sets ri(,^) (plain), r2(^) (dashed), and r3(^) (dotted) in the 
one-cut case. We have that ri(^) = [— a(^),a(^)], T2{(,) = {—ioo,—ij{^)] U 
[i7(e),ioo),andr3(C)=M. 



5.3 Proof of Theorem [531 

5.3.1 Branch points 

As a first step to the proof of Theorem 15.31 we calculate the branch points 
for the algebraic equation si{w;^) = z. 

Lemma 5.4. Let ^ > ^cr- Define n(^),v(^) > such that 



Then, the branch points are ±a{S,),±ij{(,) where 

aiO = s,{u{C),0 = 2u{0-^>0 
-iliC) = si{iv{0,0 = i (2^(0 - f^) , (5.18) 

and a(0'7(C) > 0. Moreover, we can rewrite the symbol as 

siiw,0 = y^ + + 3^3 ■ 
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Proof. The proof is straightforward. Note that ±n(^) and ziziv{S,) are the 
zeros of the derivative of si{w; ^) with respect to w. From (j5.16p and (j5.17p 
it can be shown that 3u{^f - v{^f > u{^f - 3v{^f > if ^ > ^cr- The 
positivity of a(^) and 7(^) fohows from these inequahties. D 

5.3.2 The restriction of si to R and iR 

Consider the algebraic equation 

X 3x-^ 



for real values of z. Figure 4(a) shows a sketch of the graph of the symbol 
si(x;^) for real values of x. The solutions to (j5.19p for z G M are real or 
come in pairs of complex conjugate numbers. By Lemma 15.41 ±u(^) is a 
double solution of the equation if z = iba(^). Then, as is clear from the 



graph in Figure 4(a) , the other two solutions are complex conjugate and 
their modulus is smaller than ti(^). 

Now consider the restriction of the symbol si to the imaginary axis 



3y^ 

We claim that si has four purely imaginary zeros: ziziyi,ziziy2, where yi > 
2/2 > 0. To see this, recall that n(^)^ — 3v{^)'^ > 0, so that (n(^) — 
\/3f (^))(n(^) + \/3i;(.^)/3) > 0. This can be rewritten as 

«K)^ - ..({f > MM). (5.20) 

Now consider the biquadratic equation 



By (j5.20p the discriminant is positive and less than (u(^)^ — w(O^)^ > 0- 
Then the claim follows. Figure [4(b)| shows the graph of the restriction of si 
to the imaginary axis. 

Consider the algebraic equation 

u(0^ -.■({)% MflW^^^ 

y 3y^ 

for z G M. For z = ±7(,^) the equation has four real solutions: the double 
solution =F^(0 and two strictly positive/negative solutions. The latter are 
simple and v{S,) lies between their moduli. 
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(a) Graph of x i— > s\(x\ 5). 



(b) Graph of J/ i-> — isi(ij/;^). 



Figure 4: Graphs of x i— t- S\{x\ ^) and y i— t- —isi{iy; S,), the restrictions of the 
symbol si to the real and imaginary axes. 



5.3.3 Auxiliary lemmas 

We start by proving two lemmas. 

Lemma 5.5. Let (^ > S,ct- Assume that Wa,Wb G C are such that Wa ^ Wb, 

\wa\ = {wbl, and si{wa;6.) = siiwb'jO = -2- Then z £RU iM. 

Proof. The complex numbers w'^ and w^ lie on a circle of radius p = \w'^\ = 

I'w'^l centered at the origin of the complex plane. We can factorize si{w;^) 

as 

{w^ + yj){w^ + yi) 



si{w;C) 



w^ 



where ziziyi,ziziy2 are the zeros of the symbol si(w;^). Thus 
dist(u;2,-yf)dist(u'2,-yi) 



|si(^«>6l 



p 



,3/2 



if k' 



P- 



Since — ?/i, — 2/2 < 0) it follows that 



-7r,7r] 



dist(pe*^, —yl) dist(/9e*^, —7/2) 



P' 



,3/2 



is an even function that is strictly decreasing as 9 increases from to vr. 
Thus, equality 



dist(/9e^^% -yf ) dist(pe^^", -yi) = d[st{pe'^\ -yj) dist{pe'^\ -yl) 



24 



with 9a, Ob G [—'^,M^ '^^^ only occur if 9b = ±0a- Then it follows from the 
assumptions of the lemma that w'^ = w"^ or w"^ = w"^. This gives rise to 
three possible cases: Wa = —Wh and Wa = i^Wb- Substituting these results 
into the algebraic equation yields 

z =si{wa;C) = si(-Wb;0 = -si{wb;0 = -z, 



z =siiwa;0 = siiwb]^) = si{wb;S,) = z, or 



z =siiwa;0 = si{-wb;0 = -siiwb;0 = -z. 
In all three cases z G M U iM. D 

Lemma 5.6. Let C > Ccr- Then rj(^) C M U iM for j = 1, 2, 3. Moreover, 

ri(Onr2(e) = r2(Onr3(O = 0- 

Proof. Fix j £ {1, 2, 3}. If z G ^jiC) ^^d si{w, C) = z has a double solution, 
then z is one of the branch points ±a(^) or zizij{S,)- So z £ ML) iR. If 
z G ^jiO ^iid •51 (w,^) = z does not have a double solution, then Wj{^) 7^ 
Wj-^.l{^) and |t(^j(C)l = l'"^i+i(OI- Then it follows from Lemma [5.51 that 
z G M U zM. This proves that r/(0 C R U zM. 

Now assume that for a certain value of z G M U iM and j G {1,2} 

z G r,-(0 n r,-+i(0. Then \wj{z;0\ = \wj+i{z;0\ = \wj+2{z;0\- We 
distinguish four cases. First consider the case z = 0. Recall that the al- 
gebraic equation si{w,^) = has four different imaginary solutions: Hyi 
and Hy2- This contradicts the assumption, z cannot be one of the branch 
points {±a{^),±i^{^)} either. Next assume that z G M \ {0,iba(^)}. From 
the proof of Lemma [531 it follows that 'Wj{z;^) = Wj+i{z;(,) = Wj+2{z;(,), 
so that two roots coincide and z is a branch point. This possibility was 
already excluded. If z G iM \ {0, ±i^{S,)} we analogously obtain Wj{z;S,) = 
—Wj+iiz;^) = Wj+2{z;S,). Then again z is a branch point. Since we 
could exclude all four possible cases, we can conclude that ri(^) n r2(^) = 

r2(6nr3(e) = 0. □ 

5.3.4 Transformed symbol 

Lemma 15.61 is used in the proof of Theorem 15.31 Another ingredient of that 
proof is the transformed symbol Si, defined as 

Si{W; = siHOW; = a(0 (w + ^^ + ^ + ^, (5.22) 
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see (|5.13p . (|4.6p . and ()4.7p . The advantage of the transformed symbol ^i 
over the symbol si is that it depends on ^ in a much easier way. This will 
significantly simplify the calculations. 

We denote the zeros oi^{W; by ±U{0 and ±iV{0, where C/(0, ^(6 > 
0. One can check that ±U{C) = ±u{C)/a{^) and ±iV{^) = ±iv{^)/a{C), 
where n(^) and v{^) are defined as in (|5.16p and ()5.17p . Moreover, Si gives 
rise to the same branch points as si does, namely 

Si{±U{0;0 = si{±u{0;0 = ±«(0 and 
Si{±iV{0;0 = si{±iv{0;0 = ThiO- (5.23) 

5.3.5 Proof of Theorem [57^ 

Proof. We will use the restriction of the symbol to the real axis to determine 



the sets Tj{^) n M. Figure 4(a) shows the typical form of the graph of this 
restriction. 

Recall the algebraic equation (|5.19p . For z E (— a(^),a(^)) we find two 
pairs of complex conjugate solutions. Therefore, {—a{^),a{^)) C ri(^) n 
r3(^). If z = ±a(^) the equation has the double solution =t«(^) and one 
pair of complex conjugate solutions with smaller modulus. We conclude that 
[— a(^),a(^)] C ri(^) n r3(^). Next, take z > a(^). The equation then has 
two real solutions and one pair of complex conjugate solutions. The largest 
real solution is denoted by xi{z), the smallest real solution by X2{z), and the 
complex conjugate solutions by X3(z) = X4{z). From Lemma 15.61 it follows 
that the equation does not admit three solutions with equal moduli. Thus, 
the situation |x2(2)| = 12:3(2)1 = |x4(z)| cannot occur if z > a{^). Because 
|a;2(ct(0)l > k3('^(0)l = k4('^(C))l and the roots of a polynomial equation 
are continuous with respect to the coefficients of the equation, we have that 
(a(0, +00) C TsiO and (a(0, +00) n ri(0 = (a(0, +00) n r2(6 = 0. We 
obtain similar results if 2; G (—00, —a{^)). At this moment we conclude 

ri(e)nM=[-a(e),a(e)], r2(OnM = 0, andIRcr3(e). (5.24) 

Let us now restrict the symbol si to the imaginary axis to determine 
the sets rj(,^) n iM.. Figure [4(b) | shows a typical graph of y 1— )• —isi{iy;£,). 
We proceed in a similar way. Consider for real values of z the algebraic 
equation ()5.2ip . For z G (0,7(C)) this equation has four real solutions with 
different moduh. Therefore, (0,^7(0) nFji^ = for j = 1,2,3. For 
z = we obtain the solutions ±yi,±y2- Since yi > y2, belongs to ri(,^) 
and r3(^), but not to r2(^). This is consistent with (j5.24p . If z = 7(C), 
the equation has the double solution — f (^), a real solution with modulus 
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less than w(^), and a real solution with modulus greater than t'(^). Thus, 
ij{^) belongs to r2(C). Next take z > j{C)- The equation then has two 
different real solutions and one pair of complex conjugate solutions. Using 
a similar continuity argument as before we obtain {i^{^), +ioo) C T2{C) and 
(n(0' +^co)nri(^) = {ij{^),+ioo)riT^{^) =0. The same procedure works 
for z < 0. Summarized this is 

ri(OniM = r3(e)niM = {o} and r2{0niR = iR\{-ij{0,ii{0)- 

(5.25) 
Combining (|5.24p and (|5.25p proves (|5.14p . 

To prove (a) and (b) the transformed symbol Si will be useful. First, 
we prove that ^ i— )• a(^) is an increasing function. Take ^ > ^cr- It follows 
from (iOSI) that 

Since C/(^) is a zero of the derivative of Si, the first term on the right-hand 
side vanishes. Taking the partial derivative with respect to ^ in (|5.22p yields 

MO _ f^^^^ ^ 3 ^ da(0 



Observe that the function S, ^ a(0 is increasing, see ()4.5p . Because U{S,) > 
we conclude that ^ i— )• a{(,) is an increasing function. It can be proved 
similarly that ^ i— )• 7(^) is an increasing function for ^ > ^cr- 

Next, we show that lim a(^) = +00. From ()4.5p it follows that a(^) ~ 

5-i>+oo 

V^ as ^ — )■ cxD. Using ()5.16p - (|5.18p we compute 

^ ^2(r2a(0 + e) + V(^Me)+0' + 12r2a(0') ~ V^- 



MO 
Therefore, ^ 1— )• a(^) is unbounded. The limit lim 7(^) = +00 can be 

proved analogously. 

Our final task is to calculate the limits of Oi{^) and 7(1^) as ^ — )• 0+ for 
t > t'^. In the limit ^ = the transformed symbol is 

Si{W;0) = ^ + ^„ (5.26) 
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because lim a(£) = 0. Its derivative 



r^ 



F(W^;0) = -Ti72-3TT7T 



has only two zeros, denoted by 



±iViO) = ±z\l^. 



It follows that 



/3-J-2 

lim V(0 = V(0) = \ and lim U{0 = +00. 

C^o \ t c^o 



Substituting these results into (I5.26P yields 



t r2 



lima(0 = lim5i(C/(C);e)= 1™ — + ttt^ = and 



t r2 



S'^^^^ = frnS,{-^V{0■,0 = -F(0) + HF(0)F = '"^ 



see dM]) for y*. D 

6 Asymptotic analysis in two-cut case 

6.1 Doubling the recurrence relation 

In Section Owe introduced and analyzed the symbol si{w,^) in the one-cut 
case. 

Here we want to do something similar for the two-cut case. Note, how- 
ever, that the recurrence coefficients b^n and c^^n in (|4.2p do not have limits 
as k/n — )• ^ E (0,Ccr)- Instead, there is two-periodic limiting behavior given 
by (|4.10p - (j4.13p . This is a fundamental difference with the one-cut case and, 
therefore, the construction from the previous section does not apply to the 
two-cut case. 

We analyze the two-cut case by doubling the recurrence relation (14. 2p . 
This yields a new recurrence relation in which the coefficients have limits. 
Indeed, we obtain 



x' 



'Pk,n{x) = Pk+2,n{x) + -4fc,„Pfc,„(x) 

+ Bk,nPk~2,ni^) + Ck,nPk-4.,n{^) + Dk,nPk-Q,n{^), (6.1) 
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where 

^k,n — C-kn \ Cfc-|_i ft + 0/j n^A;— l,n) 
^k,n "k,n^k—l,n ~r Cf^nOf^—^i^, 
^k,n — Cfc^nCfc— 3,n' 

The hmits of these coefficients as /c,n — )• oo such that k/n — )■ ^ G (0,Ccr) 
exist and are denoted by 

A{0= hm Ak,n = bo{0 + bi{0, (6.2) 

fc/n— >g 

5(0= hm i?fc,n = co(0 + ci(0 + ^o(0^i(0, (6-3) 

fc/n— s-g 

C(0= hm Cfc,„ = 6o(e)ci(e) + co(e)6i(0, (6.4) 

fe/n— i>g 

D(0= hm I?fc,„ = co(Oci(e), (6.5) 

fc/n— >g 

see also (liJO]l - ()ilB . 

In analogy with ()5.13p we dehne the symbol 

We use the subscript 2 to remind us that we are in the two-cut case. The 
hat refers to the fact that the recurrence relation was doubled to obtain this 
symbol. Also the quantities that are associated with the symbol (|6.6p will 
be equipped with a hat. Thus we use Wj{z;^) to denote the solutions of 
's2{w,^) = z with the usual ordering 

\wi{z;^)\ > \w2{z;^)\ > \w3{z;0\ > 1^4(2;; 01- 

Furthermore, we have '^jiC) &iid Jr, for j = 1, 2, 3. 
It is remarkable that S2 has the factorization 



S2{w; = ^""^..^^ {w^ - tr^w + t^w + r^O- (6.7) 



This follows from (J6.2p - (j6.5p and the explicit expressions for 60(0' ^i(?)> 
co(0 and ci(0 from Theorem 14.21 Note that w = —S, is always a double 
zero of (j6.7p . Thus is always a branch point and in fact an endpoint of 
one of the sets Tj{^), as will follow from the analysis in the next subsection. 
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6.2 Analysis of the symbol S2 in the two-cut case 

The remainder of this section is devoted to the proof of the following theo- 
rem, which is the two-cut case version of Theorem 15.31 



Theorem 6.1. Fix t < t'^ and < ^ < .^cr- Then we have that 

f i(C) C M+, f2(e) C R-, and f3(0 C M+. 

More precisely, there exist a(^) > /3(^) > 0, 7(^) < 0, and 5{S^) > such 
that 

r2(0 = (-00,7(0], (6-8) 

r3(0 = [?(e),+oo). 

In addition, we have for every fixed t < t"^ 

(a) ^ I— )• a(^) is strictly increasing for < ^ < .^cr with 

lim S(^) = T {t —t), and lim a{i) = lim a{^) ; 

(b) /3(^) = if and only if t < — r^ and —tr"^ < ^ < ^cr- Otherwise 
^ I— )• /3(0 is positive and strictly decreasing with 

-2/2 



\im m=T'{T'-t), 



.2 ^ i ^ _2 



lim /3(0 = 0, for -T^ <t < r' 

hm ^(e)=0, /ort<-r2; 

(c) 7(0 = if and only if t < and ^ < — tr^. Otherwise ^ 1— >■ 7(0 is 
negative and strictly decreasing with 

Af3 

lim 7(0 = -;^^ = -{y*f, forO<t< t\ 

lim 7(^) = 0, for t < 0, 

lim 7(^) = - lim 7(0^; 
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^ = -*^'\\e 




Figure 5: Regions and subregions of the t^-phase diagram. The critical 
semi-parabola separates the one-cut case region Ci from the two-cut case 
region C2. The critical ray divides C2 into the subregions C2a, C2b and C2c- 



(d) 5{^) = if and only if -r^ < t < r^ and -tr"^ < i < icv Otherwise 
^ I— )• 5(^) is positive and strictly decreasing with 



lim (5(0 = 
lim ?(0=0, 

hm m = 0, 



*\2 



(X*) 



/or t < 0, 

fort< -r^ 

/or - r^ < t < 0. 



Theorem 16.11 indicates a refinement of the t^phase diagram. We divide 
the two-cut case region C2 into three subregions C2a, C2b and C2c, depending 
on whether /3(0, 7(0 and (5(0 are zero, or not. The regions are separated 
by the critical ray 

,2 



C 



-tT' 



t<0, 
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which is tangent to the critical semi-parabola. The three regions are 



C2a 
C2b 
C2c 



t < 0, < C < -ir^ 
-t"^ <t <t'^, max(-tr^, 0) < ^ < ^c 



see Figure [5j Then, according to Theorem 16.11 we have the following for 

i < ^CT, 

. if (t,0 e C2a then ^(e) = 0, 7(6 < and 5(0 > 0; 

• if [t, ^) e C2b then /3(^) > 0, 7(0 = and 6(0 > 0; 

• if (t, e C2c then ^(^) > 0, 7(0 < and 6{^) = 0. 

The proof of Theorem 16.11 is in the following subsection. The approach 
is similar to the one used throughout the previous section, but there are 
certain complications. 

6.3 Proof of Theorem 16.11 

6.3.1 The zeros of S2 and graphs 

The symbol S2 has a double zero in — ^, see (j6.7p . The two remaining zeros 
are also negative. We order them such that xi < X2 < 0. The location of 
(t, ^) in the phase diagram determines the way the zeros are ordered, see 
Figure El 

Lemma 6.2. The zeros of's2{w;S,) are ordered as follows: 

-^ < Xi < X2, if{t,0 6 C2a, 

xi < -^ < X2, if{t,0 e C2b, 

Xi < X2 < -C, ifit,0 e C2c, 

xi = — ^ < X2, if i = —tr and t < — r , 

xi < X2 = —£,, if i = —tr and — r < t < 0, 

Xl=X2 = -e, if{t,0 = {-T^T^). 

Proof. The proof of this lemma is straightforward. D 

The three Figures [6l [71 and [H show sketches of the graph of S2 for (t,^) 
belonging to C2a, C2b, and C2C, respectively. The graphs have the following 
properties: 
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Figure 6: The graph of x i— t- S2(x;^) for (t,^) G Q 



2a- 



• S2 has four negative zeros, counted with multipHcity. — ^ is a double 
zero. 

• S2 attains a local minimum a{^) > in a point Wq > 0. 

• S2 attains three local extrema /3(.^) > 0, 7(^) < 0, and 5{S,) > in the 
respective points w^ < W2 < w^ < 0. /3(^) and d{S,) are local maxima. 
7(^) is a local minimum. Only the extremum attained at the double 
zero — ^ is zero. 

• If (t,C) G C2a, then wl = -C, so that ^(^ = 0, 7(C) < 0, and 5(0 > 0. 

• If (t,^) E C2fe, then w^ = -^, so that /3(^) > 0, 7(^) = 0, and 6(0 > 0. 

• If (t,^) G C2c, then w| = -^, so that /3(0 > 0, 7(^) < 0, and ?(0 = 0. 

6.3.2 Auxiliary lemmas 

The following two lemmas are the analogues of Lemmas 15.51 and 15.61 



Lemma 6.3. Let ^ £ (0,^cr)- Assume that WajWh G C are such that Wa 7^ 

Wb, \wa\ = \'Wh\, and 'S2{wa]i) = 's2{wb;Cj = z. Then, 2; G M. 
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Figure 7: The graph of x i— )■ 52(2;; for (t,^) G C2b- The figure shows a 
situation where /?(^) < (5(^) < S(^). It is also possible that /3(,^) > (5(.^) or 



Proof. The complex numbers itJ^ and Wf, lie on a circle of radius p = \wa\ = 

\wi,\ centered at the origin of the complex plane. We can factorize 's2{w,^) 

as 

_ {w + ^f{w-xi){w-X2) 
S2{w;0 = 

where xi,X2 < 0. Thus, 



w^ 



\'S2iw,0\ 



dist(tf , — C)^ dist(w, xi) dist(w, X2) 



if \w\ = p. 



Since xi,X2, — C < 0, it follows that 

dist(/9e*^, —S,)'^ dist(/9e*^, xi) dist(/9e*^, X2) 



[-7r,7r] ^M: 6lh^ 



p- 



is an even function that is strictly decreasing as 9 increases from to vr. 
Thus the equality for 6 = 6a, and 6 = 6b can only occur if O^ = ±6a. Then, 
it follows from the assumptions of the lemma that Wb = w^ and 



Z = S2{wb;0 = S2{Wa,S.) = {s2{Wa,0 = Z- 



Therefore, z G M. 



D 
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Figure 8: The graph of a; i— )■ S2(x;^) for (t,^) G C2c- 



Lemma 6.4. Zei ^ G (0,Ccr)- T'/ien, rj(^) C M for j = 1,2,3. Moreo?;er, 

Ti (e) n f 2(0 = r2(e)nf 3(0 = 0. 



Proof. The proof is analogous to that of Lemma 15.6 



6.3.3 Transformed symbol 

For convenience, let us again introduce a transformed symbol 52 



w 



w^ 



w^ 



see (|6.7p . Note that 52 depends on ^ in a simple way. The branch points of 
the transformed symbol 5*2 coincide with the branch points of the symbol 
S2. To see this, define W*{^) = w*{^)/^ for j = 0, 1, 2, 3. Then, 



and 



d 
dW 



52(M^*(O;O = 0, fori = 0,1, 2, 3, 



«(e) = 52(^o(0;0, m) = S2{wnc);0, 

7(0 = 52(H^2*(0;0, ^(0 = S2{W^iC);0- 



(6.10) 
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6.3.4 Proof of Theorem \6A\ 

Proof. We prove the first part of the theorem only for the case that (t, ^) G 
C2a- The other cases can be treated similarly. Figure [6] shows the graph of 
the symbol S2 for the case (t,^) G C2a- 

For z £ (7(^),0), the equation 'S2{x;^) = z has four different negative 
solutions. Therefore, (7(0,0) n fj(0 = 0, j = 1,2,3. If z = 7(0 the 
equation has the double solution w\^ a negative solution with a greater 
modulus and a negative solution with a smaller modulus. We conclude that 
7(0 S r2(0. Next, take z < ^{C)- Then, the equation has two negative 
solutions and one pair of complex conjugate solutions. Using the continuity 
argument (which was also used in the proof of Theorem 15. 3p we obtain that 
the modulus of the complex conjugate solutions lies between the moduli of 
the negative solutions. Therefore, (—00,7(0] "^ ^2(0- 

Now focus on 2: > with the extra assumption that a{^) > 5{^). We 
can make a similar reasoning if a(^) < 6{S,). If z = the equation has a 
double solution — ^ and two different negative solutions with a smaller mod- 
ulus. Therefore, belongs to ri(0- If z E (0,5(0) ^^ ^^^ two different 
negative solutions and one pair of complex conjugate solutions. The con- 
tinuity argument guarantees that the modulus of the complex solutions is 
the greatest, so that [0,(5(0) '^ ^i. Proceeding in the same way, we obtain 

[?(0,a(0]cji(0nf3(0. 

For z > 3(0 the equation has two different positive solutions and one 
pair of complex conjugate solutions. By the continuity argument we obtain 
(a(0,oo) C ri(0 or (3(0, 00) C r3(0- Since ri(0 is a compact set, see 
[7], only the second inclusion holds. 

Collecting this information we obtain 

fi(0 = [0,3(0], r2(0 = (-00,7(0], and f3(0 = [^(0,oo). 

so that (|6.8p is proved under the assumption that (t,0 ^ C2a and 3(0 > 
(5(0- One can prove all other cases in a similar way. 

Next, let us study the behavior of 3(0, /3(0, 7(0, and 5(0 as ^ increases. 
Since all branch points are of the form 5'2(W,*(0;0 for some j = 0, 1,2,3, 
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see ()6.9p and ()6.10p . we are interested in the sign of 

d - dSo dW*(C) dSo 



(6.11) 

where the last equahty holds because of (|6.9p . 

Since W^{C) = wS(C)/C > 0, we obtain from (fOTT) that ^ ^ "(0 is 
a strictly increasing function. Putting ^ = — 1 in ()6.1ip . we also get that 
that = 5'2(— 1;^) remains zero for fixed t and (t,^) in a fixed subregion 
C2a, C2b, or C2c- The other two branch points can be written as ^2{w*{^);S,) 
with i«^(^) < and w*{^) ^ —S,- In terms of the transformed symbol they 

are S2iW*{C);0 with W*{0 < and W*{C) ^ -1. Equation ([GUT) then 
implies that these branch points are strictly decreasing functions of S,- 

Let us now concentrate on the behavior of the branch points as .^ — )• ^^cr— 
for fixed t, so that (t, ^) approaches the critical semi-parabola. This behavior 
is a straightforward corollary of the following claim 

Mw^;Ccr) = si{w;^,,f, weC\{0}. (6.12) 

Let us prove this claim. Using ()5.13p . ()6.6p . and ()4.14p . we compute 

y - / 2 c\ 2 /^T , ^o(Ccr) , Co(Ccr)\ /-. , fel(Ccr) , Ci(^cr) 

= lim si(w;0^- 

The calculation of the limits of the branch points if ^ — )• 0+ can be done 
in the same way as in the proof of Theorem 15.31 and is left to the reader. D 

6.4 Reformulation of results in the two-cut case 

In what follows it will be convenient to undo the doubling of the recur- 
rence and consider instead of the functions Wj{z;^), the sets rj(^), and the 
measures fr- also for ^ G (Oj^cr); 

Wjiz;0=Wjiz'',0, j = 1,2,3,4, (6.13) 
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TjiC) = {zGC\z^eTj{0}, J = 1,2,3, 
and 

d4(z) = |.|^(z2)|d.|, zer.iO, (6-14) 

dx 

dn\(x) = \x\^(x'^)dx, xeTsiO. 
dx 

Then, with these definitions 

r,(0 = {zeC\ \wjiz;^)\ = \wj+iiz;m, J = 1,2,3; 
and 

^ 2 27ri \^u;j_(z;0 ^^^(z;^)^ 

(6.15) 
One can also check that the p.- are measures on rj(^) with total masses 

/ d/zf = 1, / d4 = g, and / d/u| = -. 
By Theorem 16.11 we have 

ri(o = [-«(0,-/3(e)]u[/3(e),a(e)], 

^2{0 = {-ico,-i-i{i)] U [i7(^),+ioo), 

r3(0 = (-oo,-<5(e)]u[,5(0,+oo). 

with 



«(0 = \/S(0, /3(0 = \/^' 7(e) = A/^7(0, and 5(0 = \/^. 
Theorem 16.11 also shows that for every fixed t < r^ 
(a) i^^ a{i) is strictly increasing for < ^ < ^cr with 



lim a(0 = ts/t'^ — t, and lim a{C) = lini 0^(0, 
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(b) 13{^) = if and only if t < -r^ and -tr'^ < C < Ccr- Otherwise 
^ I— 7- /3(^) is positive and strictly decreasing with 



lim /3(e) = T\/r2 - t, 
lim f3(n = 0, for - r^ < i < r 



2 < f - -2 

2 



lim /3(e) = 0, for t < -r , 

^-l>-tr2- 

(c) 7(e) = if and only if t < and < ^ < -tr"^- Otherwise ^ ^-^ 7(e) is 
positive and strictly increasing with 

2*3/2 

lim 7(e) = 0, for t < 0, 

lim 7(e) = lim 7(e) 

(d) (5(e) = if and only if -r'^ < t < t"^ and -tT2 < e < ecr- Otherwise 
e I— )■ (5(e) is positive and strictly decreasing with 

lim (5(e) = ^ — '^- = X*, for t < 0, 

lim (5(e) = 0, for t < -t^ , 

lim (5(e) = 0, for - r2 < t < 0. 

Figure [9] shows sketches of the sets rj(e), j = 1,2,3, for (i,e) belonging 
to C2a, C26, and Cic- 

By Theorem 15. II we know that the measures (/2^, /ig, /^l) are characterized 
by a vector equilibrium problem. The transformed measures (//^,//2,M3) 
from (j6.14p then also are characterized by a vector equilibrium problem, as 
stated in the following theorem. 

Theorem 6.5. (a) The vector of measures (/x^,/^!,/^!) given in ()6.14p is 
the unique minimizer for the functional 

Eo{pi,P2,P3) = I{Pl) -I{Pl,P2) +I{P2) - I{P2,P3)+I{P3), 

among all vectors {pi,P2,P3) of positive measures with finite logarith- 
m,ic energy I{pj) < 00, satisfying 
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r2(0 



r2(0 



ri(0 



ri(6 



r3(0 



TsCe) 



(a) 



(b) 



r2(0 



ri(0 



1^^^^^. . .L . . I ■ ■ ■ ■ ■ I . 



r3(0 



(c) 

Figure 9: The sets ri(^) (plain), r2(^) (dashed), and r3(^) (dotted) for 
{t,0 e C2a (a), {t,0 G ^26 (b), and (t,C) G Csc (c). 



(i) supp(pj) C rj(^), /orj = 1,2,3, and 
(ii) pi(ri(0) = 1, /02(r2(e)) = 2/3, and /93(r3(e)) = 1/3. 

(b) T/ie measures /Up/igj/Ug satisfy for some constant £^ 



i^-2U''^^{z) + U''kz) 



C/^i(z)-2C/'^2(z) + C/^3(2) 



[/^2(^)_2[7A'3(2) 



log 
log- 
log 



W2{z;0 
W2{z;0 



W3{z;0 
W3{z;0 



W4{z;S,) 



Proof It follows from ([O]) and (fOi]) that 



(6.16) 
(6.17) 
(6.18) 
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Applying this to the variational conditions (|5.5p ~ ()5.7p in the context of the 
symbol ?2 and using ()6.13p establishes ()6.16p - (j6.18p . D 

7 Proof of Theorem [22] 

We now give the proof of Theorem 12.21 which is based on Theorem 15.21 For 
fixed t, we define the measure vi as an average of the measures fi\ 



Jo 



yi = / lA d^- (7-1) 

Jo 

This measure will be the limiting normalized zero distribution of the poly- 
nomials pn^n as n — ;■ oo. Note that the measure ^\ is defined by (|5.4p in the 
one-cut case (i.e., ^ > ^cr), and by (|6.14p in the two-cut case. 

7.1 Proof of Theorem [Ml for t > r^ 

We first prove Theorem 12.21 for t > t"^ . This is the simplest case, since we 
deal with the one-cut case for every ^ > and we can apply Theorem 15.21 
almost immediately. 

Proof. Let t > t'^. Conditions (a) and (d) of Theorem 15.21 follow from 
Lemma 14.11 and Theorem 14.21 The interlacing condition (b) follows from 
Theorem 12. II (a), and condition (e) is contained in Theorem 15.31 So in order 
to be able to apply Theorem 15.21 it remains to establish condition (c). This 
will be done in the following lemma. 

Having verified all conditions, Theorem 15.21 can be applied and Theorem 
follows for t > r^. D 



To complete the preceding proof we still have to establish the follow- 
ing lemma. It states that the recurrence coefficients b^^n and Ck^n remain 
bounded in case k/n is bounded. 

Lemma 7.1. Let t € M. Then the two sets of recurrence coefficients {bk^n I 
k + 1 < n} and {ck^n \ k + 1 < n} are bounded. 

Proof. By ()4.3p and ()4.4p it is sufficient to proof that the set 

{ak,n I A; + 1 < n} 

is bounded. From the asymptotics of the orthogonal polynomials qn,n, see 
[H [5] , it follows that there exists M > such that the zeros of the diagonal 
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polynomials qn,n, n = 1,2,3, .. ., belong to [—M,M]. Then, by interlacing, 
see Theorem 12 -H the zeros of all q^n with k <n belong to [—M,M]. Note 
that the zeros of Qk^n coincide with the eigenvalues of the tridiagonal Jacobi 
matrix 



•Jk,n 



( y/ani 



v 



\ 



-y/«fc-l,i 
A/«fc-l,n J 



It is well-known that the eigenvalues of a real symmetric matrix X interlace 
with the eigenvalues of the matrix obtained from X by deleting the first row 
and column. Applying this k — 2 times, we find that that the eigenvalues of 







are in [—M,M], which implies that a^-i^n ^ M'^ for every k < n. This 
proves the lemma. D 



7.2 Proof of Theorem [M] for t < r^ 

The situation for t < r^ is more complicated. For small values of ^, namely 
^ < ^cr, we are in the two-cut case, while for ^ > ^cr we deal with the one-cut 
case. 

The two-cut case was handled in Section [6] by doubling the recurrence 
relation. This led to the symbol s'2 and related notions. It will be convenient 
to double the recurrence relation also in the one-cut case. The doubled 
recurrence relation is (16. ip . and so for ^ > ^cd we define 



hm Ak^n 

k/n-^S, 

hm Bk.n 

k/n-^S, 

lim Ck.n 

k/n^r^ 

hm Dk^n 



MO 
B{0 
c{0 

D{0 



k/n—^S, 



26(6, 


(7.2) 


2c(0 + 6(0', 


(7.3) 


nocio, 


(7.4) 


c(e)^ 


(7.5) 



and 






w 



w^ 



w 
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The subscript 1 is added to remind ourselves that we are in the one-cut case. 
The hat refers to the fact that the recurrence relation was doubled to obtain 
this symbol. 

As before, we consider the algebraic equation si(w; i) = z for complex z 
and denote its solutions by Wj{z\^), j = 1,2,3,4, ordered such that 

\wi{z;C)\ > \w2{z;C)\ > \ws{z;^)\ > \w4{z;^)\. 

Furthermore, we define for ^ > ^cr 

fiiO = {zeC\\wi{z;a = \w2{^;m, 

and 

ZlTl \Wl-[Z) Wl-^-[Z) J 

for z £ ri(^). The following lemma establishes the link between these 
notions and the corresponding ones related to the original symbol si. 

Lemma 7.2. Let ^ > £,cv Then, z G ri(^) if and only if z^ G ri(^). 
Moreover, 



d//f(x) = |x|^(x2)dx, x£Tx{i). (7.6) 



t 
dx 

Proof. The proof of this lemma is straightforward. We omit the details. D 

We now prove Theorem 12.21 if f < r^ . 

Proof of Theorem \KE for t < t'^. Let t < t"^. 

For even k = 21, we have that pk^n is an even polynomial and we can 
write 

Pk,n{x) = n^nix^), (7.7) 

where ri^n is a monic polynomial of degree I. Rewriting the doubled recur- 
rence relation 16.11 in terms of the polynomials r^ „ yields 

xri^n{x) =n+l,n{x)+Ai^nn,n{x)+Bi^riri-l,n{x)+Ci^nn-2,n{x)+Dl,nn~^,n{x), 

where 

Al^n = Ml^ri, Bi^n = -B2Z,„, Cl^n = C2l,n, and Di^n = D2l,n- 

Therefore, the polynomials ri^n satisfy condition (a) of Theorem 15.21 
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We have for ^ > 

lim Ak,n = A{i), lim Bk,n = B{^), 

lim Ck,n = C{0, lim A,„ = ^(0, 

k/n^^ k/n-^^ 

where A(^), S(^), C(0, and L»(0 are defined by dZ^D-dZS]) if ^ > ^cr and 
by dejD-dES]) if ^ < Ccr- It follows that 

lim Ai^n = A{20, lim ^;,„ = 5(2^, 

lim Q,„ = C(20, lim A,n = i?(20- 

This establishes the condition (d) of Theorem 15.21 Condition (c) follows 
from Lemma 1 7. II and (b) from Theorem 12. iT b). For condition (e) we need 
that ri(2^) C M. If 2^ > ^cr, this is guaranteed by Lemma [7.2l and Theorem 
15.31 If 2^ G (0,Ccr); this follows from Theorem 16.11 Applying Theorem 15.21 
we obtain 

hm v{ri,n) = 2 / /I^« d^ = / /I^ d^. 
Then by ([LZl), ([Oi]) . and ([LS]), we obtain 



lirn^ 2^(Pn,n) = lA d?> = '^i 

pvpn " U 



n even 



by the definition of vi. In a similar way one finds the same limit vi for the 
subsequence of odd n, completing the proof of Theorem 12.21 Alternatively, 
this follows from the interlacing of zeros. D 

8 Proof of Theorem [2731 

8.1 Averaging the vector equilibrium problems 

In this section we prove that the limiting zero distribution vi can be char- 
acterized by a vector equilibrium problem for three measures with external 
fields acting on the first and the third measure, and a constraint acting 
on the second measure, see Theorem 18.11 Recall that vi is given by (|7.ip . 
Similarly we define the measures V2 and v^ by 



V. 



3 



1 



^i]di, j = 2,3. (8.1) 
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It follows from Theorem 15.31 and Theorem 16.11 that the supports of fx, as 
a function of S, are increasing if j = 1,3 and decreasing if j = 2. Hence, 

rri(l)cM, fori = l, 

suppi^j) = [j supp(^f) = <^ r2(0+) C iR, for j = 2, 

0<5<i [r3(l)cM, fori = 3. 

We also know that z^i(M) = 1, z^2(*) = 2/3, us{R) = 1/3. 

The vector of measures {1^1,1^2,1^3) is characterized by a vector equilib- 
rium problem. 



Theorem 8.1. Define 



Vi{x) = 


ir- 


wi{x;() 
W2{x;^) 


d^ + 


poo 

/ ^° 


a = 


J 4de, 








Mx) = 


If'- 


W3{x;0 
W4{x;C) 


dC + 


poo 

/ lo 




dC, X G 



d^, X G 



(8.2) 
(8.3) 
(8.4) 



Then, {1^1,1^2,1^3) is the unique vector of measures minimizing the energy 
functional 

E{pi,p2,p3) = I{pi) - I{pi,p2) + I{p2) - I{p2,p3) + I{pz) 

+ Vi{x)dpi{x)+ V3{x)dp3{x), 



among all measures pi, p2 and p3 satisfying 

(a) pi is supported on M and pi{R) = 1, 

(b) p2 is supported on iR and p2{iR) = 2/3, 

(c) p3 is supported on R and P3{R) = 1/3, 

(d) p2 satisfies the constraint P2 <^- 

Proof. We have already shown that {1^1,1^2,1^3) satisfies the conditions (a), 
(b) and (c). Also (d) holds, because of the definitions ()8.ip and (j8.3p . Since 
T2{i) is a set that decreases as ^ increases, we have 

supp(fT-z^2) = r2(l). 
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Since the energy functional is strictly convex it is sufficient to show that 
the Euler-Lagrange variational conditions associated with the minimization 
problem are satisfied, namely 



2U''^{x)-U''^x) + Vi{x) 



= i for X G supp(z^i), 
> i for X G M \ supp(i^i), 



for some i, 

and 

-U''^{x) + 2U'"'{x) + V3{x) 



= for z G supp((T — 1^2), 
< for z G iM \ supp(cT — 1^2) 



= for X G supp(z^3), 
> for X G M \ supp(f3). 



(8.5) 



(8.6) 



(8.7) 



We establish ()8.5p - ()8.7p by integrating the Euler Lagrange variational 
conditions dSS]) - ([^ 



£«-2C/^i(z) + C/^2(z) = log 



W^^iz) - 2[/^2(z) + U''3(z) = log 



m'^z) - 2U''^z) = log 



wi{z;S,) 



W2{z;i) 
W2{z;i) 



wsiz;^) 
W3{z;0 



if C > ^cr and (fGTTBI) - (l6TB 



£^ - 2U''^z) + W^z) = -log 



W4{z;^) 
wi{z;C) 



U^'^{z)-2U^'^z) + U^'^z) 



C/''2(z)-2C/^3(2) 



log 
log 



W2{z;i) 

W2{z\S,) 



W3{z;i) 
wz{z]C) 



W4{z;i) 



(8.9) 
(8.10) 

(8.11) 
(8.12) 
(8.13) 



for ^ G (0, ^cr) with respect to ^. 



Integrating (j8.1ip from to min(,^cril) and (j8.8p from min(^crj 1) to 1 
yields 



i-2U''^{z)+U''^{z) 



min{^cr,l) 



log 



wi{z;i) 



W2{z;C) 



de+ / log 

min(^cr,l) 



Wi{z]^) 



W2{z;C) 



de, 
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for some constant ^ G M. 

Let X G M. Since |wi(x;^)| > |it;2(3;;^)| for every ^ > 0, we can extend 
the integration to infinity and obtain an inequality 

i-2U'''{x) + U''^{x)<Vi{x), (8.14) 

since Vi is given by (|8.2p . Equality holds in (j8.14p if and only if \wi{x; ^)| = 
1^2(3:^; 0\ fo'^ every ^ > 1. That is, if and only if x G ri(^) for every ^ > 1, 
which means 

x€ f^ri(0 = ri(i) = supp(i.i). 

The first equality holds since the sets Ti(^) are increasing as (, increases. 
This proves ()8.5p . 

The proof of (j8.7p is similar. 

Integrating (|8.12p from to min(,^cr) 1) and (|8.9p from min(^cri 1) to 1 
yields 



-^ ^min(5cr,l) 



9 / ^°S 

^ ./o 



W2iz;0 



dC+ / log 

'min(5cr,l) 



W2iz-C) 



W3{z;^) 



de. 



W3{z;0 
Since |t(;2(2;;OI ^ l'"^3(-2^;C)l ^^ every ^ > it follows that 

U''^{z)-2U''^z) + U''''{z)>0, zeC. (8.15) 

Equality holds in (|8.15p if and only if |it;2(2;;^)| = |u'3(z;^)| for every ^ G 
(0, 1). That is, if and only if 

zG fl r2(0 = r2(i) = supp(a-i.2). 

0<5<1 

The first equality holds since the sets r2(^) are decreasing as ^ increases. 
This proves (|8.6p . 

This completes the proof of Theorem 18.11 D 

8.2 Auxiliary lemmas 

From Theorem 18.11 we know that (1^1,1^25^3) is the minimizer for a vector 
equilibrium problem with external fields and an upper constraint. To com- 
plete the proof of Theorem 12. 3l we evaluate the external fields Vi and V3 and 
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the measure a and show that they are equal to the external fields Vi (up to 
a constant) and V3 and the constraint a that appear in Theorem 12.31 

Although the calculations involved may not look too elegant, we think 
it is remarkable that they can be performed at all. We start by defining and 
investigating two functions. We define for j = 1, 2, 3, 4, 

Fj{z,0 = z-T^^P--w,{z;0, e>ecr, zee, (8.16) 

where Wj{z;S,) is defined as in subsection 15.21 If t < r^ we also introduce 
for j = 1,2, 3, 4, 

G,{z,0= /i^, ^ 0<e<ecr, zee, (8.17) 

where Wj{z]^) is defined by (j6.13p . We will establish three lemmas concern- 
ing these functions. The first lemma states that Fj and Gj are antideriva- 
tives of the integrands in (j8.2p and (j8.4p . This explains our interest in these 
functions. The second lemma states that Fj and Gj continuously connect 
to each other at the boundary point ^cr of their domains of definition. The 
third lemma gives the limiting behavior of Fj(z,£,) and Gj{z,^) as ^ — t- 0. 

Lemma 8.2. Let ' denote the partial derivative with respect to z. 

(a) For every t G M and z € C the equality 

dFn ■w'Az;^ 

^iz,0 = ^j^, OCcn (8.18) 



holds. Here, Wj{z;^) is defined as in subsection \5.^ 
(b) For every t < t"^ and z G C the equality 

dGi , , 1 w'Jz;0 

holds. Here, Wj{z;^) is defined by (I6.13|) . 
Proof, (a) Let ^ > ^cr, z £ C and define Wj{z; ^) for j = 1, 2, 3, 4 as 

W,iz;C) = ^^^. (8.20) 
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Since si{wj{z,^)) = z we find by ()5.13p and the explicit expressions (j4.6 
and (jiTT]) for 6(0 and c(0 that 



t t2 



(8.21) 
see (I02|l for the definition of Si{W;0- % iKW\i and (Ig:^ . we can write 



r2 



i^:,-(^,e) = ^ - ||r(7^ - «(OH^.-(^;e), (8.22) 



so that 

^2 



Taking the derivative of (j8.2ip with respect to ^, we ehniinate a'{S,) from 
(j8.23p to obtain after some calculations 

= _^(,.e).^ I 



where the second equality follows from (j4.5p . Calculating partial derivatives 
of (I8.2ip with respect to z and ^ yields 

i i = 1 i H i = 

dW dz ' aw^ dC d^ 

from which we obtain 

dWj _ dWj dSi _ dWj _,,,,Wf + 3 



d^ dz d^ dz "^^^^ W. 



3 



where the last equality follows from (|5.22p . Now, combine this with (j8.24p 
to obtain 

dF, ^ 2M 

Recalling ()8.20p completes the proof of part (a) . 

(b) The proof of part (b) follows along the same lines. Let < ^ < ^cr 
and z€C. Define Wj{z; S.) for j = 1, 2, 3, 4 by 

Wj{z; 



W,iz;0 = ^^^^^^, (8.25) 
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so that Wj{z'^;^) = ^i^i^, see ([03|l . Then by (fSlTD we can write 

g.(^,g) = ^, ". , (8.26) 






(^'0 = -^KjOtt:^ ;^- (8-27) 



so that 



Recall the definition (j6.9p of the transformed symbol ^2 and note that 
by the definition (fOSi) 

Wj{z;0 W,iz;0' W,iz;0' 

(8.28) 

see (I6.9p . Calculating partial derivatives of (|8.28|) with respect to z and S, 

yields 

9^2 51?^ _ dS2 dWj dS2 _ 

Therefore, 

sty,- _ dWj dS2 _ dWj (Wj + 1)2 
dS, dz d(, dz w. 

where the last equality follows from (16. 9[) . If we combine this with (I8.27P we 
obtain 

This completes the proof of part (b) because of (18.250 . D 

Lemma 8.3. Assume t < r^. Then for z G C and j = 1, 2, 3, 4 the equality 
of the limits 

hm G,(z,0= lim F,(z,e) (8.29) 

holds. Here, Fj and Gj are defined as in (|8.16p and (|8.17p . 
Proof. We claim that 

lim VF,(z;e)^= lim l^,(22;0. 
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To prove the claim, observe that, as a corohary of ()6.12p . 
hm Wi(z:£) = hm Wjiz :£). 



Then, from (I4.5p it follows that 



T^-t 



hm a(0 = ^— , (8.30) 



SO that also 

.2 ,\ -2 






u5j(z2;^) 



which proves the claim. 

Note also that ([830]l and ([8:2T]l yield 

Given this, the proof of ()8.29p comes down to a calculation. Using ()8.22p 
and (I8.3ip . we can rewrite the right-hand side of (I8.29P as 

lim FAz.F) = lim ( • ; -\ -^ I = lim 



where also the second equality follows from ()8.3ip . Lemma 18.31 then follows 
from the claim and equation (j8.26p . D 

We also need the limiting behavior of the Fj and Gj functions as ^ — )• 0+. 
Here we make a connection with the function uji that is defined in subsection 
12.31 and the functions UJ2 and uj-^ that are defined on the interval [— x*,x*] 
in case t < 0. 

Lemma 8.4. (a) Ift>T'^, then 

lim Fi(x,^) = 0, for x G M, (8.32) 

lim F2{x,^) = X — TUJi{x), /or X G M. (8.33) 
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(b) Ift<T^, then 

lim Gi(x,0 = 0, for x G M, (8.34) 

lim G2ix,S,) = x-TUJi{x), forxGR. (8.35) 

If t < 0, then also for j = 3, 4, 

lim Gj{x,^) = X — TUJj^i{x), for X £ [—x*,x*]. (8.36) 

(c) Ift>T^ and z G (— ic», —iy*\ U [zy*, +ioo), t/ien 

lim (F2(z+,0-i^2(^-,e)) = 2TRea;i(z), (8.37) 

where F2{z±,^) = lim F2(2=f/i,^). 
/H-0+ 

(d) If t < T^ and z € (—ioo,—iy*]U [iy*,+ioo), then 

lim (G2(^+,0-G2(^-,6) = 2rRewi(z), (8.38) 

where G2{z±,^) = lim ^2(2; =F ^lO- 
Proof, (a) Let t > t"^. For ^ > 0, we define 

w,(z;C) = p-r, j = 0,l,2,3. (8.39) 

Since si(u;j(2;;^);^) = z, j = 1,2,3,4, see ()5.13p . we obtain that ()8.39p are 
the four solutions of the equation 

a(0 (— +3w) +tw + w^ = rz, (8.40) 

ordered such that 

\u}o{z;C)\ < \uJi{z;^)\ < \uj2{z;C)\ < l^3(^;;OI- 

As ^ — )• 0+, we have that a(^) — )■ 0. Then, (|8.40p has one solution that 
tends to as well, hence 

lim uJo{z;C) = 0, 
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while the other solutions tend to the solutions of the cubic equation 



a; +tu! = Tz, 



that we already encountered in (j2.2p . 

In terms of (j8.39p . Fj{z,^) can be rewritten as 

F,{z,0 = z-TUj^i{z;0 r^, j = 1,2,3,4. (8.41) 

Moreover, using (j8.40p we get 

T ' 

from which (|8.32p follows by letting ^ — ^ 0. 

Let X G M with x 7^ 0. Then by Theorem 15.31 (a) we have x ^ ri(^) 
for sufficiently small ^. Thus w\{x\^ and W2{x\^ are real, while ^113(2;;^) 



and W/^(x\^ are non-real and complex conjugate, see Figure 4(a) Hence, 
by (|8.39p . ijj\(x\ ^) is real for small ^ and, therefore, converges as ^ — )■ to a 
real solution of ()2.2p . In the present situation there is only one real solution, 
which was previously defined as ijJ\(x). Thus wi(a;;^) -^ ^\{x) as ^ — )• 0+ 
for X G M with x 7^ 0, but in view of continuity it also holds for x = 0. Then, 
(I8.33P is obtained by taking the limit ^ — )• in (I8.4ip with j = 1 and noting 
that a(0 -^ 0. 

(b) We prove part (b) only for t < and x G [0, x*]. The proof for the 
other cases follows from similar considerations. 
For t < 0, we recall that 

S2{wj{z;0;0 = z^ j = 1,2,3,4, (8.42) 

see Section 16.11 and (|6.13p . By Figure [5] we are in the case C2b for small 
enough ^. See Figure [7] for the graph of S2 in this case. 

Let X G [0, X*). Then x^ < 5{^) for sufficiently small ^ by part (d) of 
Theorem 16.11 It then follows that W3{x;S,) and W4(x;^) are real. We now 
distinguish two cases depending on whether x^ is smaller than t'^{t'^ — t) or 
not. 

• If x^ < T^ (r^ — t) then it follows from part (b) of Theorem 16.11 that 
x^ < l^iO fo^ small ^. Then by (|8.42p and Figure [7] we have 

wi{x;^) < W2{x;C) < -C < wsix;^) < W4{x;C) < 0. 

By the definition ()8.17p this implies 

G2(x; e) < Gi(x; < < X < G^ix; < Gs{x; 0- (8.43) 
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• If x^ > t'^{t'^ — t) then by part (a) of Theorem 16.11 we have that 
x"^ > a(^) for small enough S,- In this case in Figure [7] the local 
minimum a{(,) at Wq is smaller than the local maximum S{^) at w^. 
We then get by (f8^ and Figure [7] that 

-C < ws{x;C) < W4{x;^) < < W2{x;C) < wi{x;^). 

Thus by JKWf 

< Gi{x; < G2{x; < x < G^ix; < G^{x; 0- (8.44) 

Next, we get from ([8:i2|) and (JHTT]) that the Gj{x; , ^), j = 1, . . . , 4, are 
the four solutions of 

G^ - 3xG^ + (3x2 ^ ^^2 ^ ^^^2 ^ (_^3 ^ ^4^ _ ^^2^ _ 2^^)(- + ^^2 ^ Q_ 

Letting ^ — )• 0, we see that one of the Gj's tends to while the other three 
tend to the solutions of 

G^ - 3xG2 + {3x^ + tT^)G -x^ + t\t^ - t)x = 0. (8.45) 

In view of (|8.43|) and (|8.44p it must be Gi(x;^) that tends to and so we 
proved ([83i]l . 

Now let ujj{x]^) = ^(x — Gj+i(x;^)), so that 

G,{x-i) = x-Tujj.i{x-i), 3 = 2,3 A- (8.46) 

It then follows from ()8.43p and ()8.44p that in both cases 

a;i(x; C) > > uj^{x; i) > 0J2{x; i). (8.47) 

As ^ — ;■ 0, we have that Gj(x;^), j = 2,3,4 tend to the three solutions 
of ([8^ . Then from (|8:i6]) we get that ujj{x;^), j = 1,2,3, tend to the 
solutions of 

uj^ + tuj = Tx. (8.48) 

In view of (j8.47p and the earlier definitions of ujj{x), j = 1,2,3, it is then 
easy to check that 

lim ujj{x;^) = ujj{x). 

This proves (f835|) and (l8:36D because of (IHIiel) . 

(c) Let t > t"^ and z G {—ico,—iy*) U {iy*,ioo). Recall that uJi{z) is 
defined as the solution of uj"^ + tu = tz with positive real part. As in 
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the proof of part (a) we use again the definition (j8.39p for u}j{z;S^). These 
functions are defined with possible cuts on parts of the real and imaginary 
axis. We prove that 

lim tJi_(z;^) =uji{z). (8.49) 

where a;i^_(z; ^) denotes the limiting value on the imaginary axis taken from 
the right half plane. 

To prove (j8.49p . we note that, by Theorem IS.Sf b). z G ^2{C) for suffi- 
ciently small ^. Then, ■W2,+{z;^) and W2,-{z;S,) have the same imaginary 
part, but opposite real part. Note that W2{x;S^) is positive for sufficiently 
large x > and that W2{-;C) does not take purely imaginary values for z in 
the right half plane because it solves the algebraic equation (j5.19p . It then 
follows from the continuity of u'2(-;0 that W2,-{z;^) has positive real part. 
Then, by (j8.39p . also uji-{z,^) has positive real part for small enough ^. 
This proves (18:491) . 

Since ■W2,+ {z;S,) and W2-{z;^) have opposite real part, we know that 

F2{z+,0-F2iz-,0 = -2ReF2{z-,0 = 2rRe fwi,-(^;0 ™^^^ 



Then, (fSrsTP follows by letting ^ ^ and using ([8^ and a(^) -^ 0. 

(d) Part (d) can be proved in a similar way. We do not give details. D 

8.3 Proof of Theorem [273] 

Theorem 12.31 follows immediately from Theorems 12.21 and 18.11 and the fol- 
lowing result that connects Vi, V3 and a as defined in Theorem 18. H with Vi, 
V3 and a that are defined in subsection 12. 3[ 

Theorem 8.5. There is a constant C, depending on t and r such that 

Vi (x) = Vi (x) + C, X G M. (8.50) 

Furthermore we have 

Vsix) = Vs{x), X G M, (8.51) 

and 

a = a. (8.52) 

Proof of (|8.50p . We distinguish two cases. 
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Case 1: t > r^. In this case ()8.2p takes the simple form 



Vi{x) 



log 



W2ix;0 



de. 



Fix X G M. As ri(^) is an unboundedly increasing set, there exists ^*{x) > 
such that X G Tifc) if ^ > C(.x) and x ri(0 if C < C{x). If x G ri(e), 
we have that \wi{x;^)\ = \'W2{x]^)\, so that the upper bound in the integral 
can be replaced by ^*(x). The derivative of Vi can be written as 

f<'(^Uw[ix;() w'^ix;0' 



Vl{x) 
Part (a) of Lemma 



lim 



d^. 



Vl{x)= lim {F2{x,Ci 
6^0 



wi{x;^) W2{x;S.), 
now yields 

) - Fi(x,6)) + (Fi(x,r(x)) - F2{x,e{x))) ■ (8.53) 



By definition of S,*{x), x is on the boundary of ri(^*(2;)). Therefore, x is 
one of the branch points ibQ(^*(x)). The algebraic equation si{w, S,*{x)) = x 
then has the double solution 'Wi{x; S,* (x)) = 'W2{x;S,*{x)), so that the last 
terms of (I8.53P vanish, see (j8.16p . The first terms can be handled by (I8.32J) 
and ([8331) . We obtain 

V{{x) = X — TUJi{x). 

Using (I2.2p . we can integrate this equation with respect to x 

Vi{x) = ^- f {3uji{xf + tuji{x)) dL^i(x) = ^-^LOl{x)^-^tLOi{xf + C, 

where C is a constant of integration, which proves ()8.50p in view of ()2.5p . 



Case 2: t < r^. Vi is given by 



v^.W^ 



log 



wi{x;S,) 



W2{x;0 



dC + 



log 



wi{x;C) 



W2{x\i) 



de 



Fix X G M. The set ri(^) is unboundedly increasing if ^ increases. Therefore, 
there exists C{x) > such that x G ri(^) if ^ > C{x) and x ri(^) if 

Assume that ^*(x) > ^cr- The other possibility, C*(x) < ^cr, is simpler 
and will be left to the reader. We obtain for the derivative of Vi 



VUx) = — lim lim 
2 6-^o+C2^Scr- 



C2 



Ci 



w'i{x;C) W2{x;^) 



dC 



wi{x;(,) W2{x;S,) 

«*(^) /u;;(x;0 w'^ix;^) 
53 \wiix]S,) W2{X\C) 



lim 



de 
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By Lemma 18.21 we obtain 

y/(x)= lim (G2(x,6)-G'i(x,6)) 

+ lim lim (Gi(x,6)-i^i(x,e3) + i^2(x,6)-G2(x,6)) 

+Fiix,e{x))-F2ix,e{x)). 

As in the previous case, it follows from the definition of ^*(x) that Fi{x, S,*{x)) 
P2ix,^*{x)). The limit at ^cr vanishes as a result of Lemma [831 So what 
remains is the limit for ^i — )• 0+, which by (j8.34p and (|8.35p . yields 

Vi{x) = X — TUJi{x). 

This leads to (|8.50p in the same way as in the other case. D 

Proof of (f83T]l . From Theorem O it fohows that r3(^) = M if ^ > ^cr- 
Then, for every x S R, we have that |tt;3(x;^)| = |i(;4(x;^)|. Therefore, the 
second integral of (18. 4p vanishes. We now distinguish two cases. 



Case 1: t > 0, or t < and |x| > x*. If t > 0, then r3(^) = R for every 
^ > 0. If i < and |x| > x*, then it follows from Theorem 16. II that x G r3(^) 
for all ^ G (0,^cr)- So in both situations, we have that |tf3(x;^)| = |t(;4(x;^)| 
for every x G M and ^ G (0,Ccr)- Then, (18. 4j) ensures us that V3(x) = 0. Als 
V'3(x) = in this case, by ([ISD and ^Tfh . 



Case 2: t < and |x| < x*. If |x| < x*, there exists ^*(x) < ^cr such 
that X G r3(^) if 4*(x) < ^ and x r3(^) if .^ < Ci^), because r3(^) is an 
increasing set. Then we obtain 

^3^"^-2^^oy,, U3(x;0 mix;0) ^- 
Applying part (b) of Lemma [8]2] yields 

V^{x) = lim (G4(x,6) - G3(x,ei)) + (G3(x,r(x)) - G4(x,r(x)). 

The last two of terms cancel each other because x = ib(5(^*(x)) and G3(x, C*(x)) 
G4(x,^*(x)) by definition of ^*. The limit for ^i — )• 0+ is calculated using 
()8.36p and it follows that 

^3(2;) = t(w2(x) - ^3(3:)). 
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Integrating with respect to x yields 

V;{x) = \u2{xf + \tu2{xf - \u,{xf - \tu:,{xf + C, 

where C is a constant of integration. Substituting x = itx* one can check 
that C = 0. This proves (18311 in view of (HZ]). D 

Proof of (f832]l . By ([83]) we have 



supp(cr) = (J supp(4) = IJ r2(0 

and this is either iM in case t < 0, or (—zoo, —iy*] U [iy* , ioo) in case t > 0. 
This coincides with the support of a. 

Let z = iy with y > y*. The sets r2(i^) are decreasing as ^ increases, 
and there is a '^*(2:) > such that z £ r2(^) if and only if ^ < ^*iz). Then 
by (j8.3p we have 

da{z) _ f^" dfi^iz) 



d. 7„ d. <*«' (^^"> 

since there is no contribution to the integral for ^ > ^cr- 

The form of — j^ depends on whether we are in the one-cut or two-cut 
case. Let us assume that i < r^ and ^*{z) > S,ct so that both cases appear 
in (I8.54p . We will not give details about the other cases, which are simpler. 

Since ^*{z) > ^cr we split up the integral (j8.54p and we use (jS.lSp and 
(|6.15p to obtain 

da{z) 1 f<" (w'2_{z]0 ^+(^;0\ ^, 

— dt, 



dz 4:m Jo \w2-{z;^) W2+{z;^) 



27i"i y^,, \w2-iz;^) W2+{z;^)^ 



Applying Lemma 18.21 yields 

i^ = -i, hm (G2(z+,6)-G'2(z-,ei,)) 
dz 2tti Ci-i'0+ 

+ lim lim [G2iz-,C2)-F2iz,-C3)+F2{z+,C3)-G2{z+,C2)) 

+ F2{z-,e{z))-F2{z+,e{z)). (8.55) 

By definition of C*{z), z is on the boundary of T2{^*{z)). Therefore, z 
is the branch point i'y{^*{z)). The algebraic equation si(it;;^) = z then has 
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the double solution W2{z; ^* {z)) = 'W3{z;^*{z)), so that the last two terms of 
()8.55p cancel each other. Also the limit at ^cr vanishes as a result of Lemma 
18.31 The limit as ^i — )■ 0+ is calculated using (|8.38p which gives 

— - — = — Recoil z). 
az m 

By symmetry, the same formula is valid for z = —iy with y > y*. This 
proves the equality ()8.52p in view of the definition ()2.1ip of a. D 
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